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INJECTION AND ACCELERATION OF ELECTRONS AT A STRONG SHOCK: 
RADIO AND X-RAY STUDY OF YOUNG SUPERNOVA 201 1DH 
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ABSTRACT 

In this paper, we develop a model for the radio and X-ray emissions from Type lib Supernova (SN 
lib) 2011dh in the first 100 days after the explosion, and investigate a spectrum of relativistic electrons 
accelerated at a strong shock wave. The widely-accepted theory of the particle acceleration, so-called 
diffusive shock acceleration (DSA) or Fermi mechanism, requires seed electrons with modest energy 
with 7^1 — 100, and little is known about this pre- acceleration mechanism: We derive the energy 
distribution of relativistic electrons in this prc-accclcratcd energy regime. We find that the efficiency of 
the electron acceleration must be low, i.e., e < 10~ 2 as compared to the conventionally assumed value 
of e c ~ 0.1. Furthermore, independently from the low value of e c , we find that the X-ray luminosity 
cannot be attributed to any emission mechanisms suggested so far as long as these electrons follow the 
conventionally- assumed single power-law distribution. A consistent view between the radio and X-ray 
can only be obtained if the pre-acceleration injection spectrum peaks at 7 ~ 20 — 30 and then only a 
fraction of these electrons eventually experience the DSA-like acceleration toward the higher energy 
- then the radio and X-ray properties are explained through the synchrotron and inverse Compton 
mechanisms, respectively. Our findings support the idea that the pre-acceleration of the electrons is 
coupled with the generation/amplification of the magnetic field. 

Subject headings: Acceleration of particles - radiation mechanism: non-thermal - shock waves - 
supernovae: general - supernovae: individual: SN 2011dh 



1. INTRODUCTION 

Identifying the site of the particle acceleration, and 
the mechanism responsible for it, has been an intensively 
studied field of particle astrophysics. The most promis- 
ing mechanisms require the existence of a strong shock 
wave, exemplified by the diffusive shock acceleration 
(DSA) mechanism (or Fermi mechanism) where the par- 
ticles acquire energy through repeated collis ions between 
up- and down-streams of the shock w ave (|Fermil 1X9491 : 
iBlandford fc Ostrikerlll978l : iBelil fl978h . Supe rnova rem- 
nants (SNRs) are believed to be t he origin of a major part 
of cosmic rays up to ~ 10 15 eV (TBaade fe ZwickvlH934t 
iBamba et all 120031 : lUchivama et all 120071 ). T heacceler- 
ation of hadrons there is still controv ersial (|Ellison et all 
120001 120071: iMorlino fe Capriolill2012i ). and studying the 
relativistic electron population helps clarify the issue 
since the essence of the acceleration must be the same, 
and understanding the electron population will enable us 
to subtract the electron contribution from the observed 
high energy emissions. 

There is one key open issue in this picture - how the 
electrons are 'pre-accelerated'. For the DSA mechanism 
works effectively, a particle must have an enough kinetic 
energy exceeding a thermal energy behind the shock, in 
order to cross the shock wave and experience the re- 
peated collisions. The nature of this pre-acceleration 
mechanism should be seen in the energy distribution 
of non-thermal electrons in the MeV range, but so far 
this energy range is unexplored observationally in con- 
ventional astrophysical acceleration sites like SNRs. 

A young core-collapse supernova (CC SN) before enter- 
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ing the supernova remnant phase offers a unique oppor- 
tunity here. The SN ejecta running into the circumstella r 
material (CSM) launches a shock wave (|Chevalierlll982T) . 
At the shock wave a fraction of thermal particles are ac- 
celerated to the relativistic speed, and a magnetic field is 
generated/amplified. Relativistic electrons gyro around 
the magnetic fie ld producing synchrotron emission s ob- 
served in radio (IWeiler et al.lll989t IChevalierlll998l ). X- 
rays are frequently detected from nearby CC SNe when 
the obser vations are performe d with sufficient sensitiv- 
ities (see llmmler et all l2002t and references therein). 
For the X-ray emission mechanism, the emission from 
the thermal electrons behind the shock, the synchrotron 
emission from the relativistic electrons, and the inverse 
Compton (IC) scattering by the relativistic electrons re- 
processing the optical photons fro m the SN ejecta to the 
X-ray energy, have been proposed (jChevalier fe Franssonl 
2006b). The physical condition is similar to SNRs, but 
there is one important difference in terms of the obser- 
vational potential. The young SNe are characterized by 
the high magnetic field content at the shock front (due 
to the high CSM density). In SNRs, the typical (am- 
plified) magne tic fi eld strength (B) is ~ lOOyu Gauss 
ijBamba et ail 120031 : lUchivama et al1l2007[ ). Thus, at 1 
GHz typical of radio observations, electrons with the 
Lorentz factor 7 > 3000 are responsible for such an emis- 
sion. On the other hand, for a typical magnetic field 
B ~ 1 Gauss in young CC SNe ([Chevalier fc Franssonl 
2006b|), radio emission at ~ 1 GHz is produced by elec- 
trons with 7 ~ 100. 

A detailed model of the electron acceleration 
and the non-thermal emissions requires an inten- 
sive set of observations in the radio and X-ray 
regime. It has been done so far for type lib SN 
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(SN li b) 1993 J (iFransson. Lundavist. fc Chevalieil 
119961 : IFransson fc Bjornssonl I1998D . but the 
X-rays were likely dominated by the ther- 
mal emission in this case (|Suzuki fc Nomoto 
19951 : IFransson. Lundavist. fc Chevalie: 



1996; 



llmmlerTlV schcnba ch. fcWand 



2001 



Nvmark. Chandra. FranssonI l2009lh There is also a 



combined model for radio and X - ray p roperties of SN 
Ic 2002ap (jBiornsson fc FranssonI 120041 ) . but a quality 
of the observational data is not as great as that of SN 
1993J to be conclusive. Recently, SN lib 201 ldh discov- 
ered soon after th e explosion in M51 (|Griga et al.l 1201 It : 
lArcavi et al.l 1201 If ) provides a grea t data s et from radio 
throu gh X-rays (ISoderberg et al.l l2012t IKrauss et al.l 
l2012tlBietenholz et alj|2012t) . 

SN 2011dh was discovered in M51 (~ 8.4 Mpc) soon 
after the exp losion. It was spe ctroscopically classified 
as Type lib ([Arcavi et al.|[201ll ). which initially shows 
strong Hydrogen lines but latter develops strong He lines 
(jFilippenko et al.lll993l: lFilippenko|[l997l ). A progenitor 
star for Type lib SNe (SNe) has been suggested to be a 
massive star which has lost most of its hydrogen enve- 
lope during its evolution, either by a strong wind or by 
a transfer to a bina ry companion (|Nomoto et al.l IT993I : 
IWooslev et al.l Il99~4l ). The progenitors of SNe lib are 
suggested to be diverse in the size of the hydrogen en- 
velope, spanning from the 'extended' progenitors of the 
red supergiant dimension ('ellb' sub class, including SN 
1993J) to the 'compa ct' progenitors of the Wolf-R ayet 
star ('cllb' sub class) (|Chevalier fc So dcrberg 20l(j). In- 
vestigating high resolution images of SN 201 ldh, a pos- 
sible progenitor was reported as a n yellow supergiant 
(YSG) with the rad ius - 3OOi?0 (|Maund et al.l I20TTI: 
iVan Dvk et al.ll2011[) . However, it has not been clari- 
fied yet if the detected star is the progenitor or a bi- 
nary companion, or even an unrelated star. Indeed, SN 
2011dh has been proposed to belong to the cllb category 
based on the properties in the optical emission in the 
first fe w days (lArcavi et al.ll2011l ) and in the radio prop- 
erties (|Soderberg et al.ll2012D~ However, the theoretical 
interpretation of the early optical and radio emissions 
connecting to the size of the progenitor has not been 
completed yet, and a YSG progenitor may still be consis- 
tent with those observational constraints ([Bersten et al.l 

i2on . 

In the optical wavelength, despite the difference in the 
initial phase lasting for a few days after the explosion, 
the subseque nt evolution turned out to be very sim ilar 
to SN 1993 J (|Arcavi et al.ll201lt iMaund et al.ll201lD . In 
this phase, the energy produced by the radioactive de- 
cay chain 56 Ni — s> 56 Co — > 56 Fe is reprocessed into the 
optical wavelength. The similarity in the light curve 
and spectra to SN 1993J suggests that the properties 
of the explosion of SN 2011dh, i.e., the energy and the 
ejecta mass, should be very similar to those of SN 1993J 
(|Bersten et al.l l2012f ). despite the possible differences in 
the hydrogen envelope structure and the CSM environ- 
ment. 

Thanks to its close distance and the very early dis- 
covery, intensive follow up observations have been per- 
formed not only in optical but also in radio and X-rays. 
In this paper, we analyze the radio and X-ray properties 
of SN 201 ldh to constrain both the electron accelera- 
tion mechanism and the progenitor environment. The 



paper is organized as follow: In §2, we formulate the 
synchrotron radio emissions from the non-thermal elec- 
tron populations and apply the model to SN 201 ldh. In 
§3, we investigate the origin of the X-ray emission from 
SN 2011dh. The paper is closed in §4 with conclusions 
and discussion. 

2. RADIO EMISSIONS 
2.1. Synchrotron Emission Model Description 
We follow the formalisms developed for the radio emis- 
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sions f r om the interacting SNe by IFransson fc Biornsson 
(fl998l) :lBio7 nsson fc FranssonI (120041) (see also IChevalieij 
19981: IChevalier et al.l l2006at iChevalier fc Fransson 



2006b). An approximate expression for the optically 



thin synchrotron emission luminosity vL v is given as 



vL v ~ ttR Vn e y v 
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Here R and V are the position and the velocity of the 
forward shock, n c is the number density of the relativis- 
tic electrons. For our 'first' model, the energy distribu- 
tion of the relativistic electrons is assumed to follow the 
power law with the index of p extending down to •y ~ 1, 
following the standard prescription adopted in this kind 
of analysis. t syn is the synchrotron cooling time scale. 
Mother is the time scale for other energy loss processes 
which do not emit at the radio frequency (e.g., inverse 
Compton (IC) scattering or Coulomb interaction). The 
Lorentz factor of the electrons emitting at frequency v 
is 7„ - SOz^ 5 ^- - 5 (here v w = v/10 10 Hz and B is in 
gauss). 

Let us follow the equipartition argument for different 
energy channels, i.e., the standard prescription for the 
acceleration of the particles and the amplification of the 
magnetic fields. In this prescription, constant fractions 
(e , es) of the shock-generated thermal energy are stored 
in the relativistic electrons and the magnetic field, and 
we have the following expressions: 

B--l:l x 1()'V ( Y'_,j42- 8 (X, ) .nauss . (2) 
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Hereafter e e ,-i = e c /10 _1 and es-i = es/10" 
CSM density scale is given by A*, which is defined below 
and A* ~ 0.1 — 10 for a WR or YSG progenitor. These 
two are connected as 

'p-2- 
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where a = e c /^B and B is in gauss. 
The cooling time scales can be estimated as 
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Here .Lsn is the bolometric luminosity of photospheric 
photons from the SN ejecta. Coulomb interactions 
can be important for the CSM of high density, e.g., 
in SNe Hp and extended SNe lib (e.g., SN 1993 J) 
(jFransson fe Biornssonlfl998| ) with A* > 100. For SN 
201 ldh it is negligible, and we neglect the Coulomb loss 
in the following discussion. 

So far no assumption has been made on the expansion 
of the shock wave. We assume that the CSM density is 
well represented by the steady state wind solution, 

Pcsm = 5 x 10 n A„r- 2 . (8) 

The CSM density parameter A* is defined in a way so 
that A* ~ (A//lO- 5 M yr- 1 )(t> w /l,OOOkms)- 1 where M 
and v w are the mass loss rate and wind velocity, respec- 
tively. A typical WR star has the mass loss properties of 
M ~ 10" 5 Af Q yr" 1 and v w ~ 1, 000 km s~\ thus A* ~ 1. 
For a YSG case, M ~ 10" 6 M o yr^ 1 and v w - 20 km s" 1 , 
thus A* ~ 5. In each case, a variation in the value of A* 
spanning more than an order of magnitude is expected, 
thus it is difficult to discriminate these two progenitor 
scenarios solely from the value of A* . 

The density distribution of the outermost SN ejecta 
is well described by a power law as a function of the 
velocity (with the index n). The hydrodynamic interac- 
tion between the materials both following the power law 
distributions is well approximated by a self-similar solu- 
tion, and the propagation of the contact disco ntinuity is 
(|ChevalieHll98l IChevalier fc Franssonll2006bh . 

V c = 8xl0 9 £ 5 Y 3 (^) °' 3 V 012 i d ai W 1 . (9) 

The kinetic energy of the SN ejecta is in the unit of 10 51 
erg (E$i = _E/10 51 erg), and id is the time since the 
explosion in day. Here we adopt n — 10.2 which ap- 
proximates the oute r densi ty structure of SNe Ilb/Ib/I c 
i|Matzner fe McKed [1991 IChevalier fe Franssonl [2008h . 
One has to take it in mind that this expression may con- 
tain an error in estimating the dynamic evolution espe- 
cially since we take a specific value of n. Since the for- 
ward shock po sition is expecte d to be close to the contact 
discontinuity (jChevalierlll982ft . we assume V — V c and 
R = R c in the following analysis. 

Fo r SN 201 ldh, the CSM density is expected to be 
low (jSoderberg et al.l I2012D . and thus the synchrotron 
self-absorption (SSA) will be the dominant absorption 
process as in other SNe Ib /c and compact SNe lib 
(jChevalier fe Franssonll2006b1) . The SSA frequency is 

^SSA.io ~ 3 x 10~ 5 a 2 / 7 R 2 / J B 9 ^ . (10) 

Now we have the complete set of equations to describe 
the radio emissions from SNe. The input parameters are 
E 51 and M e j for the SN ejecta, A* for the CSM density, 
and e G and es for the microphysics at the shock wave. 

2.2. Hydrodynamic Model Result 

We can further reduce the number of free parameters 
using various constraints. The SN ejecta properties are 
well constrained by the optical emission around the peak 
(~ 20 days after the explosion): The ejecta properties 
of SN 201 ldh should be similar to tho se of SN 1993J 
(§1), i.e., E 5 i - 1 and M cj — 2 - 4 (jShigevam et alJ 



ll994HWooslev et alj[l99l . Indeed, a detailed model for 
the optical emiss ion of SN 2011dh su ggests E 5 i ~ 1 and 
M cj - 2 - 3M (jBersten et al.ll2012t) . We hereafter take 
E$i — 1 and M e j = 3M . For Lsn required to calculate 
the IC effic iency, we take a synt hesized bolometric light 
curve from iBersten et alJ ((2012) which best fits to the 
observed bolometric light curve. This is essentially iden- 
tical to use the observed bolometric light curve (e.g., Fig. 
2 of lArcavi et al.ll201lT ) except for the IC cooling effect 
at id < 5 where the observed data are missing. Now, 
by specifying the three parameters (A*, e c , es), we can 
compute the flux at any radio frequency at any epoch. 

Figure 1 shows one of the models that fit well the ob- 
served radio properties in multi bands. In this model, 
A* = 4, e = 6 x 10~ 3 , and e B = 5 x 10~ 2 (hereafter 
model A4). The effect of the IC cooling can be seen 
in the high frequency bands (16, 25, 36 GHz) around 
id ~ 20, where the observation shows small suppression 
from the adiabatic model curves as is also seen in the 
model curves with the IC cooling. 

In addition to these three parameters (A*, e c , e^), the 
electron distribution index p is a parameter, but we can 
derive this as p ~ 3 for the radio-emitting relativistic 
electrons (7 ~ 50 — 200), independently from the other 
model parameters. This is the sa me as was taken in the 
previous ana l ysis o f SN 2011dh (jSoderberg et al.l 120121 : 
IKrauss et al"]l2012f ). and also generally found in other 
SNe Ib/c and cllb (jChevalier fc Franssonl [2006bJ) . The 
radio spectral evolution of SN 201 ldh has been followed 
up to ~ 100 days when the cooling effect should be unim- 
portant, and the photon spectral index there is consistent 
with p = 3. With this value, the temporal flux evolution 
in each band, as predicted to be cx i~ 136 , is almost a 
perfect match to the observations as shown in Figure 1. 
Changing the value of p is not acceptable, as in this case 
either the spectral index or the temporal index becomes 
inconsistent with the observations. As we do not model 
the details of the intermediate optical depth to SSA, the 
deviation in the peak in each band is an artifact - in any 
case, the model is well constrained by the optical thick 
and thin phases simultaneously, and detailed modeling 
in the intermediate phase would not provide additional 
strong constraints. Except for this detail, the properties 
of the peaks in different bands are reproduced fairly well, 
supporting our assumption that the SSA is the dominant 
absorption process in the radio wavelength. 

2.3. Efficiency of the electron acceleration 

In our analysis, A*,e G , and es are no t mutually inde- 
pendent, unlike in the pr evious studies (jSoderberg et alJ 
I2012t IKrauss et al.|[2012t) . This is because we construct 
the model taking into account the optical emission prop- 
erties, using the SN properties and resulting shock wave 
evolution independently from the radio study. Our model 
requires a rather small value of e e as compared to the 
value typically assumed in the analysis of the radio prop- 
erties of SNe (i.e., e e = eb = 0.1). A question is what 
to extent the degeneracy is involved in the model pa- 
rameters to produce similar light curves and spectra. 
From the SSA peak, we constrain the combination of 

A t fe fg (jChevalierl Il998f) and from the optically 
thin flux we constrain Al 64 e c eB- The latter scaling re- 
lation is obtained by substituting equations 2-5 into 
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equation 1, for the optically thin and adiabatic phase 
assuming p = 3. Thus, e c and es can separately be 
obtained for given value of A* as long as the cooling pro- 
cesses are not important in the wavelength of interest. 
In Figure 2, we show the combination of (e c ,es) which 
produces essentially identical light curves to fit the ob- 
served light curves. The required value for e c is rather 
insensitive to A» as is evident from the scaling relations 
as described above. As an illustration, in Figure 3 we 
show model A with A* = 30 with e c and ts following 
these relations (model A30). The radio light curves are 
similar to those of model A4 (the effect of increasing A* 
will be further discussed below). 

From Figure 2, there is another important information. 
Since e# cannot exceed ~ 1/3, we set a strong constraint 
as A* > 2. With this constraint, if we further assume 
that A* - 2 - 10 as expected for SN 201 ldh (either a 
WR or YSG progenitor), then we constrain that e e must 
be in the range between 5 x 10~ 3 and 0.01, es > 0.01, 
and a must be below unity (0.02 — 1). 

The difference in our model and the interpretation in 
the p revious studies (jSoderberg et alj[2012t iKrauss et al.l 
l2012f) should b e clarified. Using the s caling relation for 
the SSA peaks. [Soderberg et al.1 ()2012l ) estimated A* ~ 3 
assuming e c = es = 0.1. They also provided the esti- 
mate with e c = 0.3, ts = 0.01 and A» = 6 to explain 
the X-ray luminosity by the inverse Compton mecha- 
nism (see §3 for details on the X-ray properties). The 
similar value s were obtained in another paper by the 
same group (jKrauss et al.l [2012V In our models, these 
sets of parameters never reproduce the observed radio 
light curves. The reason for this difference is simple: we 
do not allow the dynamics of the interacting region as 
a free parameter but rather use the one consistent with 
the optical observations. Figure 4 shows the evolution 
of the shock wav e radius in our mo del and those esti- 
mated previously ()Krauss et al.ll2012l) a s well as the con- 
strain ts from the VLBI measurements (jBietenholz et al.1 
2012). It is seen that the shock radius estimated previ- 
ously is smaller than in our model. If we fit the evolution 
of radius estimated bv IKrauss et al.l (|2012l ). we need to 
increase the ejecta mass to > 10M©, which is clearly re- 
jected from the optical properties (§1). This argument 
relies on the applicability of the self-similar solution for 
the shock dynamics (equation 9) which can be depen- 
dent on the slope of the outermost ejecta. While we 
have followed t he widely accepted recipe used in many 
literatures fe.g.. lChevalier fc Fransson1l2006bl[2008l ). it is 
possible that this description may contain an error. Still, 
we believe the difference between our result and the solu- 
tion for e c = €b = 0.1 is too large to be explained by the 
possible error in the self-similar dynamic solution. Our 
result here suggests that one must be careful to convert 
the properties of the SSA peaks in radio wavelengths to 
the real physical size if one wants to do it with the ac- 
curacy better than a factor of two; if one adopts values 
for the microphysics parameters (e c , cb) a priori, it does 
not provide a check in the internal consistency between 
the dynamic evolution and the microphysics parameters. 
The conventionally assumed values, e c = £b — 0.1, do 
not provide a solution consistent both in radio and opti- 
cal according to our analysis. 

In this sense, the VLBI measurement is very important 



since it is independent from the microp hysics parameters 
(|Chevalierlll998l : iBietenholz et al.ll2012l ). and this should 
be used as an independent and strong constraint. Al- 
though the radius predicted in our model is large, it 
is still within the standard error of the VLBI measure- 
ments. Indeed, the VLBI measurements provide another 
independent constraint on the CSM density: A* > 4 if 
we take £51 = 1 and M e j = 3M Q . 

For a demonstration purpose, we show the synthetic 
radio light curves for the 'low velocity' model B in Fig- 
ure 5. In Model B, we artificially changed the ejecta 
properties so that the radio light curves are fit with 
e c = £_b = 0.1. In terms of the self-similar solution, we 
require M c j = 10.5M Q if E 5 i = 1, and A* = 0.9 in model 
B (hereafter model Bl). Note that there is no degeneracy 
in M e j and A» if we fix the other three parameters. The 
radius is smaller than in model A4, but still larger than 
the previous SSA estimate (note that we do not have to 
obtain the same parameter set, since the coupled dynam- 
ics and the emission processes provide constraints that 
are not always satisfied by the simple SSA estimate). We 
can obtain the evolution of the radius identical to that 
given by the SSA estimate, but it requires the large value 
of A* and a different set of (e c and e_e) as shown in Fig- 
ure 2. It is seen that model Bl predicts the radio light 
curves similar to model A4 as shown in Figure 5. The 
main difference is that the inverse Compton cooling is 
more important in model Bl due to the smaller radius 
and larger photon density, but the difference is still not 
large to discriminate model A4 and model Bl purely from 
the radio light curves. If we take model B for the evolu- 
tion of the shock radius, then the constraint is placed on 
the CSM density as 0.5 < A* < 30 (Fig. 2). 

3. X-RAY EMISSION 

SN 201 ldh turned o ut to be a strong X-ray emitter 
(|Soderberg et al.l 12012) (Fig. 6). O n average the flux 
evolution in the 0.3 - 8 keV range followed the depen- 
dence of oc t , but it showed a complicated evolution 
in detail. There is a hint of the initial fast decay from 
~ 3 x 10 39 erg s _1 at 4 days to ~ 10 39 erg s _1 at 6 days. 
It was then followed by almost a flat evolution, reaching 
to ~ 2 x 10 39 erg s _1 at ~ 10 days. Around 20 days 
it started a rather fast decline, from 
~ 2 x 10 38 erg s _1 in about 20 days (i.e., oc t~ 2 
steeper). 



10 39 erg s _1 to 



3.1. Difficulties in thermal and synchrotron emissions 

There are three mechanisms suggested so far for the 
X-ray emission from SNe Ib/c and compact SNe lib, 
i.e., ther mal emission, synchrotron, and inverse Comp- 
ton (e .g., IChevalier et al.l l2006al IChevalier fc Franssonl 
I2006bh . For SN 2011dh. lSoderberg et al.l (12012ft dis- 
cussed these possibilities, and preferred the inverse 
Compton scenario. For the thermal emission, the free- 
free emission from the reverse shock likely dominates in 
the 0.3 - 8 keV range, and with A* = 4 the predicted 
X-ray luminosity is ~ 4 x 10 37 erg s _1 , nearly two orders 
of magnitudes fainter than observed. The free-free X-ray 
luminosity can be comparable to the observed luminos- 
ity only if A* ~ 30. Such a high density CSM is not 
rejected from the radio properties alone (Fig. 3). The 
IC cooling is important, but it is still consistent with the 
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observed light curves (see below for further discussion on 
the IC cooling effect). Also, because of the high shock 
velocity, the free- free absorption time sc ale is estimated 
to be a few days (|Chevalier et al.ll2006al) . i.e., still negli- 
gible compared to the SSA and thus consistent with the 
observed radio properties. 

However, following arguments make this interpreta- 
tion disfavored: (1) the mass loss rate is then an or- 
der of magnitude larger than expected from the typi- 
cal mass loss rate of a WR or YSG star. (2) Related 
to this argument, assuming the microphysics parameters 
(e c and es) are shared by other SNe Ib/c, then it will 
lead to the increase of the CSM density (i.e., A*) by 
a factor of ~ 30 in other SNe Ib/c as well - although 
the high CSM may apply to SNe lib and in particular 
SN 2001 ldh with the possible YSG progenitor (which 
may require a special binary channel as a progenitor sce- 
nario), it is hard to reconcile such a high density CSM 
with the WR progenitor generally accepted for SNe Ib/c. 
(3) The X-ray light curve shows a complicated tempo- 
ral structure, while the free-free emission predicts a sim- 
ple power law decay unless the CSM structure deviates 
significantly from a smooth steady-wind solution. Such 
a variation in the CSM structure sh ould be seen as a 
modulation in the radio light curves dRvder et all 12004 
iSoderberg et al.ll2006t iWellons et al.ll2012D . but it is not 
seen in SN 201 ldh. 

The synchrotron emission is even less likely. Phe- 
nomenologically, if we extrapolate the radio spectrum to 
the X-ra y band then vL v (X-ra y) ~ vL v (radio) ~ 10 37 
erg s _1 (|Soderberg et a.1.1 l2012h . The situation is even 
worse if we consider the cooling processes. At 10 days, 
the synchrotron cooling frequency is at most < 10 13 
Hz. The Compton cooling frequency is less model de- 
pendent, and it must be below ~ 10 11 Hz even if a 
very strong B ~ 10 gauss is assumed. Thus, the in- 
verse Compton cooling is important to reduce the num- 
ber of high energy relativistic electrons emitting in the 
X-ray frequency. With this effect taken into account, 
the synchrotron X-ray luminosity must be many order 
of magnitudes smaller than observed. Although X-ray 
emission from some SNe Ib/c (typically those in a later 
phase than considered here) could be attributed to the 
synchrotron emission if the intrinsic relativistic electron 
spectrum flattens at the energy hig her than probed by ra- 
dio em issions (e.g., at 7 ~ 1000) ([Chevalier fc Franssonl 
2006b), this mechanism would not operate for the early 
X-ray emission at ~ 10 days as observed for SN 201 ldh. 
The reason is evident from the above argument about the 
cooling processes. Indeed, the Compton cooling energy 
of the electrons is rather insensitive to any underlying 
model assumption (i.e., it depends on the SN luminosity 
and the shock radius, but independent from the micro- 
physics parameters), and it is 7 ~ 150 — 200 for SN 
2011dh at 10 days. Thus, even if the spectral flatten- 
ing appears in the higher energy, because of the cooling 
the expected X-ray emission would not exceed the value 
estimated with the extrapolation of the radio spectrum. 

3.2. Inverse Compton mechanism 

The remaining, and most attractive, possibility is the 
IC scattering of the SN photospheric photons by the rel- 
ativistic electrons. We can estimate the inverse Compton 
luminosity by a formula similar to the synchrotron radio 



emission: 



nR 2 Vn e j 2 ' 



1 



(ii) 



The SN photospheric emission has the characteris- 
tic energy of ~ 1 eV, thus the electrons with 7 ~ 
A/lkeV/lcV ~ 30 up-scatter the photospheric photons 
to ~ 1 keV. On the other hand, the characteristic en- 
ergy of the electrons corresponding to the IC cooling fre- 
quency is not sensitive to the model details, and it is 
Tic ~ 150 — 200 at 10 days. Thus the electrons emitting 
the X-rays through the IC mechanism is in the adiabatic 
phase. 

In this paper we show that the IC is indeed a favored 
interpretation. Altho ugh it has been also favored by 
ISoderberg et all (|2012ft . our conclusion is based on a dif- 
ferent line of arguments than theirs. Figure 6 shows the 
X-ray light curves from Models A4 and Bl, both fit the 
radio light curves fairly well. The predicted luminosity is 
about 25 to 100 times fai nter than observe d , as is consis- 
tent with the analysis by ISoderberg et al.1 (|2012[ ) . How- 
ever, given the striking similarity of the IC prediction to 
the observed qualitative behavior as well as the difficul- 
ties in the other mechanisms, we investigate if there is 
a solution to remedy the problem and provide a unified 
explanation for the radio and X-ray properties. 

It has been suggested that changing the micro- 
scopic parameters, assuming effective electron acceler- 
ation (e c ~ 0.3) and inefficient magnetic field genera- 
tion (eg ~ 0.01) will d o to fit the X-ray data by the IC 
(jSoderberg et al.ir2012D . However, this option is rejected 
from our analysis, which constrains these parameters to 
the opposite direction (Fig. 2). A second argument is 
from a consideration on the electron number density. As- 
suming the standard single power law distribution, at 10 
days the number density of relativistic electrons assumed 
in this formalism is n rc \ ~ 6.4 x 10 5 (e c /0.1)(A*/4), while 
the total electron number density is n c ~ 1.3 x 10 5 (A*/4) 
(including thermal electrons). Therefore, the condition 
e c < 0.02 must be met, otherwise the model is self- 
inconsistent. 

Finally, even if we ignore these discrepancies and al- 
low a large value of e c , such a model is in principal not 
accommodated with the radio properties. For the IC 
scattering works effectively, the low value of B < 0.5 
Gauss is required. Then, since the energy of the electron 
suffering the significant IC cooling is model-insensitive 
and tic ~ 150 — 200 during 10 — 30 days, this inevitably 
means that the radio emission above v > 20 GHz is sig- 
nificantly suppressed by the IC cooling effect. This is not 
seen in the radio data - the observed radio spectral in- 
dex does not show the expected considerable steepening 
at the high frequency (note that there is a hint of the 
small IC cooling effect there, but this behavior is com- 
patible to that predicted by our standard model A4: Fig. 

To highlight the final point, in Figure 7 we show an 
example of the radio light curves for a model in which 
the IC cooling is effective enough to explain the X-ray 
luminosity (Fig. 6). Since model A does not allow a 
large value of e e , we start with model B and change A„ 
so that we reaches the X-ray luminosity comparable to 
the observations, while fixing the microphysics param- 
eters so that the radio emission without the IC cooling 
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fits the radio light curves (i.e., along the constraint on 
e c and es in Fig. 2). Figure 7 shows the model ra- 
dio light curves. We require A* ~ 20 (i.e., model B20). 
This model predicts the X-ray light curve through the IC 
mechanism fairly well (Fig. 6), but then the IC cooling 
effect becomes significant and alters the spectral index 
of the radio emission at > 20 GHz. It never fits to the 
observed radio properties. This is not a model-specific 
problem, but a general problem in the models with large 
value of e e since such a model inevitably requires small B 
and the resulting radio emitting electrons' energy above 
the IC cooling energy. 

3.3. What is missing? - the pre- accelerated electron 
population 

As shown in §3.2, the IC scattering cannot explain the 
X-ray luminosity in a self-consistent manner. However, 
given the failure of the synchrotron mechanism and the 
required high value of A» in the thermal scenario, and 
also the striking similarity of the IC X-ray model curve 
cvolutioin to the observed one, we further investigate 
what is required if the X-ray is indeed a consequence 
of the IC mechanism. 

All the three problems in the IC scenario found in 
§3.2 are indeed related to one single assumption in our 
(and conventional) model prescription - the single power 
law distribution of the relativistic electrons. Alterna- 
tively, we suggest that this assumption is not appro- 
priate in the low energy regime, and that the discrep- 
ancy in the predicted and observed X-ray flux does re- 
flect the deviation of the electron distribution from this 
standard/conventional formalism. The strong constraint 
from the radio properties is obtained only for the elec- 
trons with 7 ~ 50 — 200. On the other hand, the elec- 
trons responsible for the X-ray emission through the IC 
mechanism have the lower energy, about 7 ~ 30. Thus, 
a consistent view on the radio and X-ray properties of 
SN 2011dh can only be obtained if the electron number 
density at 7 ~ 30 is by about two orders of magnitudes 
larger than the extrapolation from that at 7 ~ 50 — 200. 
Adopting this distribution, the X-ray light curve is well 
reproduced (Fig. 6). The energy range of this additional 
electron population is below the energy probed by the 
radio properties, thus the resulting radio light curve is 
hardly affected and is consistent with that derived as- 
suming the single power law distribution. This relatively 
low energy electron population must be peaked around 
this energy regime (7 ~ 20 — 30) and cannot extend down 
to 7 ~ 1, otherwise it violates the electron number con- 
servation. The derived distribution is shown in Figure 
8. We suggest this population represents the electrons 
accelerated by the 'pre-acceleration' mechanism. With 
this distribution of relativistic electrons, the radio (Fig. 
1) and X-ray (Fig. 6) emissions are simultaneously ex- 
plained. 

4. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the radio and X-ray 
properties of SN lib 2011dh. From the radio light curves 
in multi-bands, we have obtained strong constraints on 
the efficiency of the electron acceleration (e c ) and the 
magnetic field generation/amplification (cb) behind the 
shock wave. 



The model developed here is only weakly constrain- 
ing the CSM environment. From the requirement that 
es < 1/3 and the shock radius through the VLBI mea- 
surement, we have a strong constraint as A* > 4. The 
fact that we see only a slight effect of the IC cooling 
in the radio light curves places a constraint that a can- 
not be much larger than unity: In the model B sequence 
A* ~ 20 is clearly rejected, while in the model A sequence 
A* ~ 30 is still acceptable. The allowed range is consis- 
tent with both the WR and YSG progenitor scenarios. If 
A* ~ 30, the X-ray luminosity can be explained by the 
thermal emission from the materials behind the (reverse) 
shock wave although the fit is rather poor to explain the 
time evolution. Although the required CSM density ex- 
ceeds the typical mass loss rate of a WR or YSG star 
by an order of magnitude, a variant of a binary inter- 
action scenario to produce the SN lib progenitor may 
(although does not have to) lead to s uch a high mass 
loss r ate in the end of the evolution (|Benvenuto et al.1 
l2012f ). The radio data alone do not reject this possibility 
for a particular case of SN 201 ldh. However, if we nat- 
urally extend our argument that the microphysics of the 
magnetic field amplification and the particle acceleration 
is similar between SN 201 ldh and other SNe Ib/c, then 
the required CSM density for other SNe Ib/c should in- 
crease by a factor of ~ 30 and never fit to the Galactic 
WR mass loss rate. Thus, we have investigated another 
possibility - how to explain the radio and X-ray prop- 
erties simultaneously for the standard mass loss rate of 
A* ~ 1 - 10. 

The IC cooling effect is important in constraining the 
properties of the microphysics and the mass loss (as men- 
tioned above for the upper limit of A*). The large values 
of A* and resulting large e c are rejected because such a 
situation will predict the significant IC cooling in the ra- 
dio wavelength. We thus obtained that e c < 0.01, much 
lower than conventionally assumed for the SN-CSM in- 
teraction. S uch a low value of e e was also derived for SN 
ellb 1993J (jFransson fc Bjornssonl I1998D. and also no t 



rejected for SNe Hp in general ( Chevalier et al.ll2006a|) . 
Thus, it can be a property shared by SNe Ib/c and 
cllb as well. Since es can as large as conventionally 
assumed (i.e., es = 0.05 in our standard model A4), the 
mass loss rate estimated thr ough the SSA scaling (e.g., 
I Chevalier fc Franssonll2006bf ) will not need a drastic re- 
vision. 

This argument on the IC cooling effect led us to con- 
clude that the X-ray properties of SN 2011dh cannot be 
accommodated by the IC mechanism, as well as the ther- 
mal and synchrotron, as long as the single power law dis- 
tribution of the relativistic electrons is assumed. Alterna- 
tively, we suggest that the single power law distribution 
is not appropriate for the energy range probed by radio 
and X-rays from SN 201 ldh (i.e., 7 < 200), and that the 
discrepancy does reflect the deviation of the electron dis- 
tribution from the conventionally assumed power law. To 
explain the radio and X-ray properties simultaneously, 
we need a population of relativistic electrons peaking 
around 7 ~ 20 — 30 in addition to the power law com- 
ponent extending to the higher energy (in which those 
with 7 ~ 50 — 200 are probed by the radio observations). 
This low-energy population indeed occupies a major part 
of the relativistic electrons both in number and energy. 
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We suggest this electron population represent the elec- 
trons accelerated by the 'pre'-acceleration mechanism, 
and the power law component represents ones further 
accelerated by the DSA-like mechanism. Although the 
argument is based on modeling radio and X-ray proper- 
ties of SN 201 ldh, we speculate that these properties may 
be shared by other SNe as well, given that SN 2011dh is 
a canonical SN in its radio properties and that the en- 
ergy range of electrons probed in the radio and X-rays in 
other SNe cllb/Ib/Ic should be similar to the case of SN 
201 ldh. Good X-ray data are generally missing in other 
SNe Ilb/Ib/Ic, and the future intensive observations in 
radio and X-ray will be valuable to study the low en- 
ergy relativistic electron population as we suggest in this 
paper. 

The typical energy of the pre-accelerated electrons we 
have found is 7 ~ 20 — 30, thus the number density of the 
relativistic electrons at 10 days is n e ~ 10 4 cm' 3 . This 
requires ~ 1 — 10% of the thermal electrons suffer from 
the pre-acceleration, rather than ~ 100% of the electrons 
accelerated to the relativistic energy required when as- 
suming the single power law electron distribution. Since 
this reduction in the number density is compensated by 
the increase of the energy content of each electron by the 
typical energy 7 ~ 20 — 30, the 'intrinsic' efficiency of the 
acceleration including this pre-acceleration population is 
similar to the 'apparent' efficiency of e c < 0.01 derived 
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by the radio modeling alone. 

Our finding provides an observational constraint on 
the nature of the pre-acceleration mechanism. The typ- 
ical energy scale, 7 ~ 20 — 30, suggests that the pre- 
acceleration of the electron takes place simultaneously 
and coupled with the amplification of the ma gnetic field 
(jFrederiksen et al.ll2003t iHededal et all 12004) . The pro- 
cess would equalize the kinetic energy of a thermal pro- 
ton a nd a (pre-)accelerated electron (jFre derikse n et al] 
I2003T) . thus 7 - (m p /m c )(V/c) 2 ~ 20 (where c is the 
speed of light, and V ~ 0.1c is the velocity of the shock 
wave). This is exactly the energy scale we need for SN 
2011dh. Then, about 1% of these pre-accelerated elec- 
trons are further accelerated through the DSA-like pro- 
cess, as evidenced by a power law behavior of the electron 
distribution in the higher energy as derived by the radio 
properties. 
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Fig. 1. — Multi frequency radio light curves of Model A4 (red solid) as compared with those of SN 2011dh. The parameters are 
(A t , e e , £b) = (4, 6 X 10 — 3 , 5 X 10~ 2 ). The synthetic light curves without the inverse Compton cooling are also shown (blue dashed). 




Fig. 2. — The values of e e and tg derived for SN 2011dh through the radio light curves, as a function of A* 
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Fig. 3. — Multi frequency radio light curves of Model A30 (red solid) as compared with those of SN 2011dh. The parameters are 
(At,,e c ,e B ) = (30, 0.01, 9.7 X 10~ 4 ) (Fig. 2). The synthetic light curves without the inverse Compton cooling are also shown (blue dashed). 




Fig. 4. — The evolut ion of the radius of the i nteraction region. Thick gray line is that e stimated with the scaling relation for the SSA 
peaks (Chevalier 1998) assuming e c = eg = 0.1( Soderbcrg et al. 2012; Krauss ct al. 2012). The VLBI measurements are given by filled 
squares) Bictcnholz ct al. 2012). Model A is shown by red filled circles (A4: A* = 4) and by red open squares (A30: A* = 30), while model 
B is by blue filled circles (Bl: A* = 0.9) and blue open squares (B20: A* = 20). 
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Fig. 5. — Multi frequency radio light curves of Model Bl (red solid) as compared with those of SN 2011dh. The parameters are 
(j4* , e e , 6b ) = (0.9,0.1,0.1). The synthetic light curves without the inverse Compton cooling are also shown (blue dashed). 
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Fig. 6. — X-ray light curves of SN 2011dh (0.1 - 8 keV) as compared with some models discussed in the text. Dashed curves are for 
the free-free emission, with A* = 4 (thin) and A* = 30 (thick). Thin lines are for the inverse Compton scattering, for model A4 (A» = 4) 
and model Bl (A* = 0.9). The inverse Compton emission in model B20 (A* = 20) is shown by the thick blue line. Finally, model A4 but 
introducing another electron population peaking at 7 ~ 25 is shown by thick red line. 
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Fig. 7. — Multi frequency radio light curves of Model B20 (red solid) with large e c to reproduce the X-ray light curve by the inverse 
Compton scattering. The parameters are (A t ,e c ,e B ) = (20,0.26,2.5 X 10~ 4 ). The synthetic light curves without the inverse Compton 
cooling are also shown (blue dashed). The effect of the inverse Compton cooling is important in the radio wavelengths. 




Fig. 8. — Constraints on the energy distribution of relativistic electrons at the shock front of SN 2011dh. The scale is for 10 days since 
the explosion. The dashed line is a standard single power law in model A4 to fit the optically thin radio synchrotron emission. Along the 
line, the red region is the constraint from the radio light curves in the optically thin phase. The red region with the label 'IC is the electron 
population required to explain the X-ray emission at ~ 1 keV by the IC mechanism. For comparison, we show the electron spectrum in 
model B with A* = 20 (B20) that explain the X-ray emission by the IC (but failed to explain the radio properties). This model B line is 
similar to the 'IC population derived based on model A, demonstrating that the argument is rather model-independent. The gray box at 
7 ~ 1 shows the thermal electron population for an appropriate range of the CSM density and the shocked radius. Finally, a schematic 
picture is shown for the distribution required for the low energy electrons to explain both the radio (by synchrotron) and the X-ray (by 
the IC), labeled as the 'injection spectrum?'. 



